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Abstract. 

Some combinatorial properties of fixed boundary rhombus random tilings with 
octagonal symmetry are studied. A geometrical analysis of their configuration 
space is given as well as a description in terms of discrete dynamical systems, 
thus generalizing previous results on the more restricted class of codimension-one 
tilings. In particular this method gives access to counting formulas, which are 
directly related to questions of entropy in these statistical systems. Methods and 
tools from the field of enumerative combinatorics are used. 

Key-words: Random tilings - Generalized partitions - Configurational entropy - 
Discrete dynamical systems - Young tableaux. 

Introduction 

The experimental discovery of quasicrystalline alloys led to extensive work on 
space tilings over the last 15 years, as it became clear that quasiperiodic, and not 
only periodic, structures could play important role in solid state physics. Indeed, 
the atomic structure of the highest quality quasicrystals has been found to follow 
closely the 3-dimensional icosahedral analogues of the celebrated pentagonal Penrose 
tilings 1^. Among the many questions that are still open in this field, the origin of 
their stability is one of the mostly highly debated. Physical explanations range from an 
electronic stabilization mechanism (refinements on the old Hume-Rothery approach) 
to an original entropic stabilization, allowed by specific phason modes which can be 
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generated in quasiperiodic tilings. Our purpose here is not to discuss the relative 
merits of the different mechanisms, but to analyze in detail the combinatorial problems 
associated with configurational entropy in random tilings. 

This paper follows a previous one Q in which the general framework was 
introduced, as well as specific results concerning codimension-one tilings. The d- 
dimensional random tilings of interest are made of rhombi {d—2) or rhombohedra 
(c?=3), or even higher dimensional analogues. These tilings are projections onto a 
d-dimensional Euclidean space of a c?-dimensional faceted membrane cut into a D- 
dimensional hypercubic lattice {D > d). The "codimension" of a tiling is the difference 
D — d, and the the tiling is said to be of type D ^ d. In reference 1^ , we discussed the 
codimension-one case for tilings with specific (fixed) boundary conditions. This allows 
us to write a one-to-one correspondence between tilings and combinatorial objects, 
called partitions. We built a geometrical description of the partition configuration 
space in terms of integral points in a high dimensional space, the entropy being 
computed from the integral volume of a specific convex polytope in that space. The 
occurrence of multiplicative and additive formulas for this volume was analyzed in 
detail, and given a simple geometrical meaning in the latter case in terms of a simplicial 
decomposition of the convex polytope. 

The aim of the present paper is the analysis of random tilings of higher 
codimension, starting with the simplest 4^2 case. Studying these cases is of 
direct importance in the context of quasicrystal physics, since all the quasiperiodic 
tilings encountered in this field are of codimension greater than one (5^2 for the 
pentagonal Penrose tiling, and 6 ^ 3 in the icosahedral case). Tilings of type 4-^2 
correspond to the so-called octagonal family, which was also observed in concrete 
alloys Q. Although they are the simplest, "octagonal" random tilings already present 
most of the difficulties which, up to now, have forbidden the derivation of exact results 
for the large class of random tilings derived from hypercubic tilings Q . Note that exact 
results exist (for the entropy) for other kinds of tilings, such as the square-triangle 
tiling ^ 0, rectangle-triangle tilings ||, ||], or large codimension tilings Note 
however that the present point of view does not apply to the two first examples since 
there exist no partition representation for such tilings. 

Our analysis for the 4^2 tilings follow from a generalization of the simple 
partition problem, valid in case of codimension one, to an iterated partition problem, 
which was proposed earlier [ ^t] , and has already led to some preliminary numerical 
results. Here we describe the intricacy of the configuration space, which is no 
longer convex "as a whole" , but remains convex by parts. We show that despite 
its complexity, some exact but partial enumerative results can be obtained, although 
we must stress that the ultimate goal - an exact formula for the entropy - was not 
obtained and seems out of reach for the moment. We nevertheless believe that the 
present analysis is an important step in at least two directions: we give a very precise 
description of the configuration space and its simplex decomposition and we point out 
several very closely related problems in combinatorics, like the enumeration of sorting 
algorithms (Appendix B). 

The paper is organized as follows. Section 1 recalls some older results and 
definitions, in particular the concept of de Bruijn lines and faceted membranes, and the 
bijection between standard partitions and codimension-one tilings. Section 2 focuses 
on higher codimension tilings, by introducing "generalized" partitions, and describing 
the particular structure that is inherited by the configuration space. Its properties in 
terms of local rearrangements of tiles (flips) are analyzed in detail. In section 3, we 
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discuss the decomposition of the configuration space into normal simplices, and we 
show the latter can be characterized thanks to a "descent theorem" . This allows us to 
compute new enumerative formulas which were inaccessible by "brute-force" methods; 
these formulas are displayed in section 4. 

Even though this paper focuses on two-dimensional tilings and more precisely on 
octagonal ones, some results can easily be generalized to higher dimensional systems. 



The state of the art in the I? — > cZ cases is briefly discussed in Appendix D 



1. Definitions and known results 



In this paper we consider 2-dimensional tilings of rhombic tiles which fill a region of 
the Euclidean space without gaps or overlaps. The standard method for generating 
such structures consists of a selection of sites and tiles in a 4-dimensional lattice 
according to certain rules, followed by a projection onto the 2-dimensional subspace 
along a generic direction. We then say that we have a 4 ^ 2 tiling problem, or an 
octagonal one, in reference to the sub-class of ideal quasiperiodic Ammann tilings 
which have octagonal symmetry. The above procedure is also known as the "cut-and- 
project" method By construction, the so-obtained rhombic tiles are the 

projections of the 2-dimensional facets of the 4-dimensional hypercubic lattice. There 
are 6 different species of tiles, two squares and four 45 degree rhombi. Figures I and 
pi] show examples. In the cut-and-project language the difference between the higher 
and the lower dimensions is called the tiling codimension. In this case it is equal to 2. 

We first recall some definitions and results which will prove to be useful 
throughout this paper. These definitions are given in a slightly more general context 
than the octagonal case. The higher dimension will be denoted by D and the lower 
one by d. 



1.1. De Bruijn grids and directed membranes 

Firstly, it should be mentioned that there exist two related classes of objects which 
can be put in one-to-one correspondence with random tilings: de Bruijn grids on the 
one hand, and directed membranes on the other hand. 

De Bruijn grids ^] are dual representations of tilings which can be useful to 
state or prove some results concerning tilings. There are a great number of publications 
dealing with these grids in the scientific literature (for example, see ^ ) , therefore 
we shall not give a complete presentation of these objects. Instead we shall give them 
an intuitive definition in the case of two-dimensional tilings. De Bruijn grids are made 
up of lines, the so-called de Bruijn lines, which are also called "worms" . These lines 
join together the middles of opposite edges of rhombic tiles. Since the tiles are rhombi, 
it is always possible to extend these lines through the tiling up to the boundary. Such 
lines are displayed in figure ^ Any tile is crossed by two lines. There is no triple 
intersection point (condition of regularity). On the other hand, there are lines which 
can never intersect, even in an infinite tiling. They join rhombus edges of the same 
orientation, as illustrated in figure |]. We say that these lines belong to the same 
family. A family is in correspondence with an edge orientation. In a D — > 2 tiling, 
there are D edge orientations and therefore D families of de Bruijn lines. 

The relevant object here is not the grid itself but the underlying intersection 
topology, which defines the tiling: a grid can be directly read on a tiling by joining 
together the middles of opposite edges, but it can afterwards be continuously deformed 
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Figure 1. A patch of octagonal (4 — ► 2) tiling. Some worms (de Bruijn lines) 
arc represented. There are 4 families of worms. We have only drawn lines out of 
one of them (grayed) . 

provided no triple point appears in the process. This grid and the tihng are said to 
be dual In the following, we will sometimes distinguish between the terms "worms" , 
which are sequences of rhombi of a tiling, and "de Bruijn lines" , which are elements 
of a grid in an abstract grid space, with no underlying tiling any longer. 

Conversely, it can be proven that, given such a grid, it is possible to build a unique 
tiling, the de Bruijn grid of which is identical to the grid under consideration [p^ |l6|. 




Two lines of two different families can but need not intersect. A grid where all 
lines of all families intersect is said to be complete. In this case, to insure the existence 
of all intersections, we impose that, "far" from the intersection region, the lines are 
perpendicular to vectors Ui , one per family. The lines of a given family are therefore 
parallel at the infinity. 



Directed faceted membranes are representations of tilings in hypercubic lattices of 
higher dimensions, which have been developed to study random tilings in parallel with 
the partition method (see below) |l9[ |ll|, H, |, |2l|, They arc the generalization 
of one-dimensional directed walks (or polymers) in hypercubic lattices. This point of 
view is closely related to the cut-and-project method. Therefore we shall only give a 
brief presentation of these membranes. The main idea is that a 13 — *■ d random tiling 
can be lifted as a d-dimensional non-flat structure embedded in a D-dimensional space. 

This structure is a continuous membrane made of d-dimensional facets of the Z 
hypercubic lattice. When this membrane is projected along the suitable direction, the 
projections of these facets are precisely the tiles the tilings are made of; its continuous 
character guarantees the absence of gaps in the so-obtained tiling. Such a membrane 
is said to be directed to emphasize the fact that its projection does not create any 
overlap. For example, figure |^ displays a 3 ^ 2 tiling, which can also be seen as a 2- 
dimensional non-flat directed membrane embedded in a cubic lattice. To get a tiling, 
this membrane must be projected along the (1,1,1) direction of the cubic lattice. 
This point of view can be generalized to arbitrary dimensions and codimensions. This 
correspondence is always one-to-one. 
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Figure 2. 3-dimonsional representation of a 3 ^ 2 tiling. 
1.2. Partitions in codimension one 

It is possible to derive from this membrane representation a coding of random tilings 
by combinatorial objects called "partitions" |2l[ |. 

This point of view is easily understood when looking at figure |^: the membrane 
can be seen as a stacking of unit cubes in 3 dimensions and an integral height (the 
number of cubes) can be assigned to each of the kl columns of this stacking, resulting 
in a X ^ array containing integers. Since the original membrane is directed, these 
numbers are decreasing in each row and in each column of this two-dimensional 
array. This latter array is called a plane partition and each integer a part. In this 
representation, the integer p is called the height of the partition. It is the upper 
bound of each part. There is a one-to-one correspondence between such partitions and 
membranes embedded in a fc x ^ x p piece of cubic lattice. There is a straightforward 
generalization of this point of view to D -t- 1 — > D membranes and ZJ-dimensional 
partitions (called hypersohd partitions), which are families of integers arranged in 
D-dimensional arrays, decreasing in each direction (for more complete details, see 
reference [||, sections 2.1 and 2.2). 

In the following section, we generalize this partition point of view to any 
codimension tilings, which enables us to build their configuration space. This general 
point of view was only briefly tackled in previous references |l|, |2^, ^ . It was developed 
and formalized in reference pl[ |. 

2. Higher codimensions tilings 

In this section, we show how it is possible to code octagonal tilings, or more generally 
D — > d tilings, as generaUzed partitions, that is families of integral variables, but living 
on structures more complex than the previous rectangular arrays. These structures 
will turn out to be the dual graphs of relevant rhombus tilings. 

2.1. Generalized partitions 

Our goal in this section is to prove that D ^ d tilings can also be coded by "generalized 
partitions on {D — I) d tilings". Let us explain what this terminology means. 
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Generally speaking, we define a partition problem as a family of K integral 
variables, denoted by Xi,X2, ■ ■ ■ ,xk, placed at the vertices of a directed graph, so 
that any two variables placed at two adjacent vertices satisfy an order relation in 
agreement with the orientation of the edge between those vertices^. The underlying 
directed graph is called the base of the partition problem. To simplify, we shall 
consider that all order relations are weak (xi > xj). The integral values are between 
and an integer p, called the height of the partition problem. A solution of this 
problem is called a partition, of height p. The integral variables Xi are called the 
parts. Figure ^ displays an example. In reference we mainly studied hypersolid 




Figure 3. A simple example of directed graph. It defines a partition problem. 
The associated partition problem has 5 variables, which are related by: xi > X2 ; 



partitions, the graph of which is equivalent to a piece a hypercubic lattice, in the 



context of codimension-one partition problems (see above, section 1.2). 

To introduce the tiling coding by partitions, we shall work in the grid 
representation. We focus here on the D — > 2 case (the presentation of the general 
D d case would require some more definitions and refinements. The interested 



reader will refer to references |21|, see also Appendix D). Let us consider a. D ^ 2 
grid. We single out a family of lines, which can be chosen as the _D-th one without 
loss of generality. It contains kjj de Bruijn lines. The D — 1 remaining families define 
a new grid. We call it a subgrid of the first one. Our goal is now to build a partition 
on this subgrid that codes the initial tiling: a part will be attached to each vertex of 
this subgrid. 

Firstly, we need to introduce the so-called interhne indices. Since they do not 
intersect, the fcu singled out lines divide the plane in fc^i + 1 domains. These domains 
are unambiguously labeled from to fc^) in the simplest way: two adjacent domains 
are labeled by two successive numbers which are increasing in the direction of ujj 



(as defined in section 1.1). Now the value of the part attached to a subgrid vertex 
is simply equal to the interline index of the domain in which this subgrid vertex lies. 
The maximum height of these parts is kjj. 

There is a more simple way of characterizing the order between these parts: since 
a de Bruijn line of the subgrid is transverse to all the lines of the _D-th family, the 
parts on this line are ordered in the same direction as the interline indices. Therefore 
we have defined a "canonical" order on every subgrid line. We say that we have 
ordered those de Bruijn lines. By convention, we chose those lines to be ordered in 
the direction of decreasing parts (we insist on this point because it is a source of 
confusion). Now, since any two adjacent vertices of the subgrid are joined by such 
a line, we have ordered any two parts. Therefore we have defined on this subgrid a 
partition problem of height kjj. To sum up, we have coded any D ^ 2 grid as a pair: 

X To begin with, we shall suppose this directed graph to be acyclic, that is to say there is n o sequence 



of inequalities such as Xi-^ > Xi^ > ■ • ■ > Xi^ > . The general case will be discussed in A-ppendix 



D 



Fixed-boundary octagonal random tilings: a combinatorial approach 



7 



a Z? — 1 — > 2 subgrid and a partition on it. 

Conversely, given such a pair, the _D — > 2 grid from which this pair comes can be 
easily re-constructed. One must add the D-th family of lines in such a way that all 
the vertices of the subgrid lie in the interline, the index of which is equal to the part 
attached to this vertex. The constraints on the parts insure that we actually obtain a 
D — > 2 de Bruijn grid. 

In conclusion, we have derived a one-to-one mapping between D ^ 2 grids and 
partitions on Z) — 1 — > 2 subgrids, the parts being suitably ordered on oriented de 
Bruijn lines. A more mathematical formulation, related to this work, can be found in 
reference 

This mapping can be translated in the tiling (or directed membrane) language: 
the generalized partitions can be defined on the suitably oriented dual graphs of the 
corresponding D — 1 ^ 2 tilings. For short, we call them "partitions on tilings". 

Figure ^ provides an example of 4 — > 2 tiling seen as a partition on a 3 — > 2 tiling. 
The 3 — > 2 tiling has been slightly deformed to anticipate the next step of the process. 
Note that parts are ordered on each de Bruijn line (or worm). Once the partition 
has been chosen, zones where parts are equal are separated by bold lines, which are 
"opened" to form worms (shaded) of width 1. This step is the manifestation in the 
tiling representation of the fourth de Bruijn family. 




Figure 4. A 4 2 tiling coded by a generalized partition on a 3 ^ 2 tiling. 
Left: the inequalities between the tiles. Right: a partition of height 2 and the 
corresponding 4 — > 2 tiling. It fills an octagon of sides 2, 2, 2, 1. 

To conclude this paragraph, we shall say that a, D d tiling problem can be 
studied as a collection of partition problems on a set of (dual graphs of) D — 1 ^ d 
tilings. Even though we have only proved this point in the 4 — > 2 case, the 
demonstration can be generalized Practically, to build a D — > d tiling, one 

can iterate a partition-on-tiling process. The first step is simply a codimension-one 
partition on a d-dimensional hypercubic array. It generates a d -I- 1 — > d tiling. The 
next steps increase D by one each. Therefore there are D — d steps. 

As compared to usual random tilings, partition-generated ones have specific 
polygonal boundary conditions. For example, the tiling in figure ^ have a hexagonal 
boundary. In the case of 4 ^ 2 tilings, the polygon is an octagon of sides fci, ^2, 
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and (see figure |^). More generally, such tilings have zoi30topa_Q boundaries. Note 
that they are dual to complete de Bruijn grids. It should also be mentioned that such 
polygonal boundaries have a strong macroscopic influence on tilings, which results in 
a lower entropy than in free or periodic-boundary systems |^ |2^, ^ . 

2.2. Configuration space 

In this section, we study the configuration space of partition-generated tilings that fill 
a given polygonal domain. 

The codimension-one case has already been studied in detail : the configuration 
space C consists of all the integral coordinate points (integral points) lying into the 
convex polytope defined by the system of inequalities related to the partition problem. 
This configuration space is embedded into an Euclidean space of dimension K, where 
K is the number of parts of the partition problem. Two points are neighbors in C (i.e. 
they are linked by an edge of the underlying hypercubic lattice) if they only differ 
by a local rearrangement of tiles which is usually called an (elementary) flip |^ (see 
figure 1). 

In this section, all the latter properties are extended to generalized 4 — > 2 
problems, in particular to partitions-on-tiling problems. In reference [ pT| , the general 
D — > d case is treated. 

Let us consider 4-^2 tilings which fill an octagonal region of sides fci, ^2, ^3 and 
ki. They are described by a class of partition problems on 3 ^ 2 tilings inscribed in 
hexagons of sides ki, k2 and k^. These tilings will be indexed by an integer a. They 
have exactly K = fcifc2 + ^1^3 + ^2^3 parts. Therefore each configuration space Ca 
related to the partition problem on a is of dimension K. Now to describe the whole 
configuration space of the tiling problem, we need to make explicit how these different 
Ca are connected to each other. 

Firstly, we need to specify how a fiip in the tiling representation is translated in 
the grid space. It is simply a 3-linc flip, as illustrated in flgure |^. 




Figure 5. An elementary flip (left) and its grid space counterpart (riglit). 

If the fourth family of lines has been singled out in the partitions-on-tiling process, 
two cases must be distinguished: 

+ Tlie zonotope generated by tlie family of d-dimensional vectors V2, ■ ■ ■ , vjj) is the set 

Z= l^a.Di, < Qi < l| . 

It is also called the Minkowski sum of vectors Vi. In two dimensions, it is a 2i5-gon. The link between 
vectors Vi and the D-dimensional representation is specified in references [pi and 
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I: either the 3-line flip does not involve any line from the fourth family It means 
that the 3 vertices involved in the flip have the same interline index and that 
this index does not evolve during the flip. These vertices are therefore coded by 
parts of same value. On the other hand, the base tiling (i.e. the tiling dual to 
the 3-family subgrid) on which the partitions are defined undergoes a flip; 
II: or the flip involves a line / of the fourth family. In this case, the 3-family subgrid 
is not modified through the flip. The same holds for the base tiling on which 
the partitions are defined. On the other hand, let us consider the only vertex S 
involved in the flip but which does not belong to I. During the flip, its interline 
index is increased by ±1. Therefore a part (and only one) of the partition problem 
varies (by ±1). 

To sum up, a flip is translated either in a base tiling flip, without any modiflcation 
of the parts (type-I flip), or in a variation of one of the parts without any modiflcation 
of the base (type-II flip) . 

Let us go back to the conflguration space C. Since C can be seen as a collection of 
spaces Ca associated with tilings a, it can be given a "discrete fiber bundleQ' structure, 
the base B of which is the configuration space of (base) tilings a. Its fibers are the 
spaces Ca- 

Practically, suppose that B is embedded in an hypercubic array of dimension 6 
and that the dimension of each fiber is K. Then the whole space can be embedded 
in a lattice Z'^ x Z^^: the first S coordinates code the base tilings a and the K last 
ones the parts on these tilings. We already know the structure of C inside a fiber: an 
edge between two vertices corresponds to a type-II flip. Now, we must establish how 
a type-I flip connects two flbers. 

A type-I flip consists of a flip in the base tiling, transforming the tiling a into a' , 
but which does not alter the values of the parts. In the flber, the K coordinates of 
the corresponding points are therefore unchanged. But in the base B, a and a' are 
coded by two points which differ by only one of their d coordinates. Thus, in C, the 
two tilings differ by only one coordinate: they are neighbors in Z*^ x Z^. 

However, we have omitted to deal with a subtlety in the previous statement: so 
far, we have proven that two fibers are connected via a piece of hypercubic lattice. 
Thus we have only proven the local hypercubic structure of the configuration space. To 
provide a complete proof, we need to exhibit an extrinsic|| set of hypercubic coordinates 
in which every configuration can be encoded and in which two neighbor tilings differ 
by a single flip. As a matter of fact, we only have to specify coordinates in fibers: the 
choice of coordinates in the base B is an irrelevant question. As it was stated above, 
a choice of coordinates is equivalent to the choice of a tile-labeling of a 3 — > 2 base 
tiling. Now, as illustrated in figure |^, such a tiling can be seen as a domino tiling on 
a triangular lattice: every tile is the union of an upward and a downward triangle. 
Therefore any labeling of upward triangles will provide a tile-labeling and therefore a 
set of coordinates in each fiber. It is now clear that with such coordinates, a type-I 
flip corresponds to a bond of the hypercubic lattice. 

Remark: we can now derive the dimension of the hypercubic lattice in which C is 
embedded: 

dc = S + K = 2kik2 + kiks + k2k3, (1) 

* We employ improperly this term. In particular all the fibers are not necessarily identical, 
(t That is independent of the fiber. 
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Figure 6. A 3 — ► 2 tiling as a domino tiling on a triangular lattice. Each tile is 
labeled like its unique upward triangle. 



since the dimension of the base is ^1^2 and the dimension of fibers is fcifc2 + fcifc3 + fc2fc3- 

These arguments in the octagonal case can be extended by induction to the general 
D — > c? problem pl[, at least as far as the local hypercubic structure is concerned 



(see Appendix D 



Note also that, following this work, this configuration space has been recently 
investigated further by M. Latapy, who detailed the nature of its structure when it is 
seen as a partially ordered set in any D ^ 2 case. Briefly speaking, the lattice 
structure in the fibers looks like that of a plane partition configuration space (it is a 
"distributive lattice" ) , whereas the structure of the whole set is a bit less rich (it is a 
lattice, but not a distributive one). 

Another interesting question concerns the connectivity of this configuration space: 
for a given boundary, is it possible to obtain a random tiling from any other one via a 
sequence of elementary flips? In the D ^ 2 case, the configuration space is connected 
for any D ^ . The present analysis provides another straightforward proof of this 
result: every fiber is connected as the configuration space of a partition problem on 
an acyclic directed graph since it is the convex union of normal simplices. Moreover, 
the base is connected for the s ame reason (i nductively), which completes the proof. 
The general case is discussed in Appendix D| . 



3. Decomposition of the configuration space into normal simplices: 
general case 

3.1. Simple descent theorem 

We first recall the results of references ^ about generalized partitions. The 
configuration space of a generalized K-pavt partition problem of height p is embedded 
in a i^-dimensional Euclidean space, the coordinates Xi of which are the parts of the 
problem. The configurations are coded by integral-coordinate points (called integral 
points), which belong to the convex polytope JFI^I defined by the intersection of the 
hypercube {0 < Xi < p) and of the cone [xi < Xj) defined by all the suitable relations 
between the parts. 

The key point is that this configuration space can be decomposed into elementary 
volumes, the so-called normal simplices. Let (ei, 62, . . . , bk) be the orthonormal basis 
of the Euclidean space which generates the Z'^ lattice. A if-dimensional simplex of 
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vertices Aq, Ai, . . . , Aj^ is said to be normal if there exists an integer s such that: 

• each Ai is an integral point, 

• AiAi+i is parallel to a vector Cfc. for any i, 

• if i j then fc^ 7^ kj , 

• IIAjAj+ill = s for any i. 

For short, we shall call such a simplex a normal simplex of side s. Its integral volume 
{i.e. the number of integral points it contains) can easily be derived: it is the binomial 
f s + K \ 

coefficient ( ^ I . But such simplices have lower-dimensional faces in common 

that contain integral points, which must not be double-counted. Therefore we must 
take into account a subtle inclusion-exclusion scheme in order to count correctly the 
number of configurations. To sum up, one must suppress j faces (j = 0, . . . , jmax) to 
some simplices in order to avoid double-counting. If aj is the number of simplices that 
lose j faces, then the number of configurations is 

(2) 

where the maximum number of suppressed faces, jmaxj depends on the partition 
problem under study. 

The so-called descent theorem ^ provides a prescription to characterize the 
coefhcients aj. To state this theorem, we need the following definition: with each 
simplex of the decomposition, we can associate the sequence (fci, . . . , kx) of indices 
appearing in the definition of the normal simplex. Then the number of descents in 
this sequence is the number of indices such that ki > ki^i. 

The descent theorem states that, if there exists a zero-descent simplex, which is 
always true up to a re-indexing of the basis vectors^, then a simplex with j descents 
loses j faces. Therefore the coefficient Oj in equation ^ is equal to the number of 
simplices with j descents. 

As a corollary, the number of normal simplices in the decomposition is equal to 
the sum of the coefficients Cj of equation ^. 

These coefficients Oj can be given a different equivalent interpretation j^, ^ : one 
builds a directed graph, denoted by T, with two extremal vertices, O and 5*0. A simplex 
of the decomposition is put in one-to-one correspondence with maximal walks in T 
{i.e. going from O to Sq). More precisely, to each vertex of the graph, it corresponds 
a configuration of height p = 1 of the partition problem. Two configurations are 
neighbors if they differ by only one part Xi, which is in the "lower" configuration 
and 1 in the "higher". Therefore the link between the two configurations can be 
indexed by i and the descent theorem can be translated in terms of these indices: the 
number of descents of a walk in the graph is defined as the number of descents of the 
sequence of indices of the bonds it follows. 

To sum up, in the graph T of any generalized partition problem, a step between 
a vertex and one of its neighbors in a maximal walk amounts to increasing one of the 
parts from to 1. Therefore a maximal walk from O to So amounts to a labeling of 
the parts, from 1 to K, which specifies in which order they are increased from to 1. 
O (resp. 5*0) is the configuration where all the parts are equal to (resp. 1). 



t At least in the case of an acyclic base graph (see A.ppcndix D ) 
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In codimension-one D + l ^ D partition problems, we have proved that the graph 
T is the configuration space of the D ^ D — I partition problem on a hypercubic array 
of sides fci, fc2, . . . , kD_i [||. 

In codimension larger than one, that is in the case of partition-on-tiling problems, 
the parts are attached to the tiles of the D — 1 ^ d problem, as in figure ^. Therefore, 
to each maximal walk in the graph T, it corresponds a labelling of the tiles, which 
characterizes in which order the parts are increased by one in the walk. 

For example, figure ^ (left) shows a tile labeling (among many others) in the 
partition problem of figure (left). In fact, the only condition on those labelings is 
that when two tiles Xi and xj are adjacent, if the order relation is Xi > Xj, the label 
associated with Xi is smaller that the label associated with Xj {xi is increased before 
Xj). In other words, these labels are ordered on de Bruijn lines. 



3.2. Decomposition in simplices 

A 4 — > 2 tiling problem is a collection of generalized partition problems on 3 — > 2 
tilings. On each such tiling, the descent theorem can be applied. Therefore, the 
counting polynomial of the 4^2 problem, which is the sum of all the individual 
polynomials on each 3-^2 tiling, can also be written 

where K is the number of tiles, independent of the 3^2 tiling, and M is the greater 
of the integers jmax involved in the collection of partition-on-tiling problems. 

In this section, we shall prove that the above result for codimension-one 
problems can be generalized: the sum of the coefficients aj of the counting 
polynomial W is equal to the number of maximal walks of a given class in the related 
configuration space of a given 3 — > 1 tiling problem. At the end of this section, we 
shall give an explicit analytic expression of this number of walks |^l| . A general result 
concerning walks in the configuration space can be derived in the general D d 
case. However, to avoid inessential complication, we shall present it in the restricted 
4^2 case, and in an informal manner. A rigorous proof in the Z3 — > 2 case is given 



in Appendix A and the general case is discussed in Appendix D 



Let us consider a partition problem on a 3 — > 2 membrane, or, equivalently, on a 
3^2 tiling. Figure (left) provides an example. 

As we have seen it at the beginning of this section, the sum of the coefficients 
aj of this (generalized) partition problem is equal to the number of labelings of the 
tiles, with integers running from 1 to K , which respect the following condition: these 
labels must be increasing on each oriented de Bruijn line. Figure || (left) displays such 
a labeling in the case of the 3^2 tiling of figure 0. 

On the other hand, a configuration of this partition problem of height p = 1 is 
characterized by the line which separates O's and I's on the tiling. Figure ^ (right) 
shows a partition of height 1 and the corresponding line. In figure ^ (right), some 
other such line configurations, associated with partitions of height 1, are displayed. 
These lines are directed walks using three kinds of elementary steps: north, west and 
north-west ones. The number of steps in each direction is determined by the side 
lengths of the hexagonal boundary. 
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Figure 7. A partition problem on a 3 ^ 2 tiling. Left: > symbols show the 
order relations between tiles. We recall that parts are canonically ordered on each 
worm. Right: a solution of height p = 1 and the line separating O's and I's. Note 
that this line has been extended up to the top-left and bottom-right corners of 
the hexagon. As a matter of fact, this line is the projection of a 3 — ► 1 tiling. 




Figure 8. Left: a walk in the configuration space of the partition problem on 
the 3 — ► 2 membrane encoded by a labeling of the 8 tiles of this membrane. This 
labeling respects the order on the de Bruijn lines. Right: the walk in the 3 — > 1 
problem configuration space associated with the left-hand-side labeling. One goes 
from a broken directed line to the next via a single flip. The extremal tilings are 
the first and the last ones. 



Now, a walk in the space of partitions of height 1 is also a walk in the space 
of these directed lines. For example, the walk of figure || (left) is encoded into the 
sequence of lines of figure ^ (right). On the other hand, we notice that these latter 
lines can be seen as projections in two dimensions of 3 ^ 1 directed membranes (i.e. 
one-dimensional walks embedded in a cubic lattice). These membranes lie on the same 
rectangular parallelepiped as the initial 3^2 membrane. A walk counted by the sum 
of the coefficients aj of the partition problem on this membrane is therefore also a 
walk in the configuration space of a given class of 3 — > 1 membranes. 

To sum up, the sum of the coefficients aj of this partition problem is equal to the 
number of ways of labeling the partition parts with some rules. Such a labeling can 
be put in correspondence with walks in a suitable configuration space of 3 — > 1 tilings. 
And the same holds for each individual partition problem on a 3 — > 2 tiling. 

Conversely, given such a walk it is always possible to reconstruct the 3^2 tiling 
it comes from, as well as the labeling on this tiling. Therefore this walk is counted by 
the sum of the coefficients aj of a partition problem on a 3 — > 2 tiling, hence by the 
sum of the coefficients aj of the initial 4 — > 2 tiling problem. 

In conclusion - and we rigorously prove this result in Appendix A - the sum of 
the coefficients aj of a 4 — > 2 tiling problem is equal to the number of maximal walks 
in the configuration space of a suitably related 3^1 tiling problem. 
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Theorem: For any problem of enumeration oj fixed boundary 4^2 tilings, the sum 
of the coef 



icients Oj of the additive counting polynomial 



M 



K 
K 



(4) 



is equal to the number of maximal walks in the configuration space of the associated 
3 — > 1 tiling problem (as defined in Appendix A). 



Generalizing codimension-one notions, by "maximal walk" , we mean a walk 
between two tilings which play a singular role in the configuration space, the so-called 
"extremal" tilings (figure |8|) . One of the properties of these extremal tilings is that 
they go in pairs and that the (Manhattan) distance between two such configurations 
is a local maximum in the co nfiguration sp ace, even in the general D — > d case pi| . 
These t ilings are des cribed in Appendix A and Appendix D (see figures Al and pi]) . 

In Appendix B, we explicitly calculate this number of walks in the Z? — > 1 



configuration spaces. In particular, if A:i{k,l,m) denotes the sum of the coefficients 
aj of the polynomial counting the number of 4 — > 2 tilings inscribed in an octagon of 
sides (fc, I, m,p), 



A^{k, I, to) = {kl + km + /to)! 



[(/-1)!M]^ 



(2/ + fc + TO-l)![2l 

(fc- 1)![21(to - l)![2l(2/ + fc - 1)!W(2/ + to- 1)!^ 
(2/-l)![2l(/ + /c-l)![2](/ + m-l)![2l ' 



(5) 



where we have used the generalized factorial functions of order 2, as defined in 

k 

reference fcl'^l — JJj!. Some properties of these functions are given in 

reference [||. 



3.3. Generalized descent theorem 

As it was done in reference Q in the codimension-one case, we shall now refine the 
previous theorem in order to characterize among the maximal walks those which are 
counted by a given coefficient Oj . 

In other words, we are looking for a descent theorem in any codimension. As 
in codimension one, we need an edge labeling in the associated 3 — > 1 configuration 
space such that Oj is the number of maximal walks which have exactly j descents with 
respect to this labeling. 

Now, we know that such an edge in the configuration space is in one-to-one 
correspondence with a facet of a 3 ^ 2 tiling (the value of which changes from to 1). 
This two-dimensional facet also belongs to a membrane attached to the rectangular 
parallelepiped P. Therefore it belongs to the piece of cubic lattice bounded by P. 

If we choose a labeling of all these two-dimensional facets bounded by P, it will 
induce a labeling of the facets of any 3 — > 2 membrane attached to P. Such a labeling 
will be used to index the variables of the partition problem on this membrane. For 
example, figure ^ displays a labeling of the facets bounded by P, and figure |l^ shows 
the induced labeling on a membrane attached to P. Then, according to the descent 
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theorem, if there exists a zero-descent walk for this labeling, then a coefficient aj^ 
of this individual partition problem counts the number of jQ-descent walks on this 
membrane, and thanks to the correspondence between facets and edges, a jp-descent 
in the whole 3 — > 1 configuration space. 

Conversely, let us consider any jo-descent maximal walk in this configuration 
space. This walk is counted by the sum of the coefficients Oj of the partition problem on 
an individual 3 — > 2 membrane. If there exists a zero-descent walk on this membrane, 
this walk is more precisely counted by the coefficient aj„ . 

In conclusion, we are looking for a labeling of the facets which induces, on any 
3 — > 2 membranes, a labeling with a zero-descent walk. Then the number of jq- 
descent walks in the 3—^1 configuration space will be equal to the coefficient a.j of 
the polynomial W{p). We propose such a labeling in the following paragraph. 

Let P be a rectangular parallelotope of sides fci x fc2 x fca embedded in a 3- 
dimensional lattice, the orthogonal basis of which is (61,62,63) (see figure^; among 
all the possible choices, the hexagonal non-flat frame of the membranes is chosen as 
follows: it contains the vertices (0, 0, ^3) and (fci, fc2, 0) of P). Each facet is coded by 



«3 




Figure 9. Labeling of tlie two-dimensional facets of a parallelepiped P of sides 
3 X 3 X 1. For any fca larger that 1, the same labeling could be handled in the 
same way. 

the coordinates of its center. We define on these coordinates an order relation -< close 
to the lexicographic order, apart from a slight difference: let /i and /2 be two facets, 
the center coordinates of which are respectively {x\,x\,x\) and {x\,x\,x\). Then 

• if ^3 > xl then -< p. 

• \i x\ — x"^ and this value is an integer, i.e. these facets are horizontal^, then 

- if < xl then -< p. 

— li x\ = X2 and if < xf then ^ p. 

• if = and this value is a half-integer, i.e. these facets are vertical, then 

f Normal to e-j. 
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- if < xl then -< p. 

— li x\ = x\ and if x\ > x\ then p -< p (only this point differs from the 
lexicographic order definition). 

Figure ^ illustrates this point. It displays a facet labeling compatible with the 
previous order relation. P has been "exploded" to enable a better reading. Let us 
emphasize the difference between horizontal and vertical layers. 

In the following figure |lO|, we have represented a 3 ^ 2 membrane attached to 
P and its facet labeling induced by the previous one. We check that this labeling 
is compatible with the order on the de Bruijn lines (which are oriented from top to 
bottom here). This result does not depend on the chosen 3 — > 2 membrane. 



A 




Figure 10. Example of labeling induced on a 3 — > 2 membrane. 



Note that the so-obtained labelings differ from those previously defined in order 
to prove the hypercubic character of the configuration space (section p.2[ ) . 

This establishes the existence of a descent theorem in the octagonal case of 
interest. It will be used in the following section to get finite tiling enumerations. 

Remark: It is natural to wonder if there exists such a result in the general D ^ d 
case. Except the facet labeling, the same arguments hold in this general context: 
if such a labeling existed, it would provide a descent theorem. However, we have 
good reasons to believe that such a labeling does not exist. Indeed its existence 
would mean that in the general case, the counting polynomial can still be written 



^(-P) ~ ( ^ ' X ^ ^ ' know examples (in the 7 — > 3 case) where 



p + K - j 



such an expression is false (see Appendix D). Therefore if there exists a general 



descent theorem, it will have a more complex statement than the previous one. 



4. Enumeration results 

^.1. Two exact formulas 



In reference 28 , Elnitsky provides two exact formulas for 4 — > 2 tiling enumeration, 
in the case where two sides of the octagon are equal to 1. 

As displayed in figure ^l], two cases must be considered, according to whether the 
two sides of length 1 are adjacent or not: 



w^xl.- E E ( )( T )( '^t' )( )^(«) 



a-ffc— r c+d— s 
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Figure 11. Two tilings of octagons of sides r,l,s,l (left) and r',s,l,l (right). 



2{r + s + l)\ (r + s + 2)! 
r! s! (r + 2)! (s + 2)! ' 



(7) 



The sketches of the proofs of these formulas are recalled in Appendix C.l 



In his paper, Elnitsky also says that no closed formula is known for the first 
double sum. We have somewhat simplified the first formula: 



(r + s + 1)! 



2(r 



^ 2fc- 1 



k=0 



r,i,s,i rl sl {2r + l){2s + 1) 
Since the last sum can be written in terms of a hypergeometric function 



k=0 



(9) 



this result is a "closed" form for the enumerating function, even though it is not 
written in terms of products or ratios of simple functions. 

In order to prove equation |^, we need the following recursion relation, due to 
Brock M: 



W, 



4^2 



r + s 
r 



(10) 



The sketch of the proof of this recursion relation is also given in Appendix C.2 . 

If the right-side member of equation ^ is written as the sum of two terms, we 
denote the first one by Ar_s and the second one by Then 



r + s 
r 



-^r.s ^r— l,s ^r,s— 1 

by simple algebraic manipulations, and 

r + s\ 1 
r ) (4r2 - l)(4s2 - 1) 



1 - 4rs 



(4r2 - l)(4s2 - 1) 



(11) 
(12) 



1 

1) ^2fc-l 



fe=0 



(2r - l)(2s - l)(r + s + 1) - r(2r + l)(2s - 1) 



k 

r ~ k 
r 



-s(2s + l)(2r-l)- 
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r + s \ 4rs — 1 f r \ f s 



r 



2 

r + s 



(4r2 - l)(4s2 -1) \ k J \ k 

4rs - 1 



r J (4r2 - l)(4s2 - 1)^ 



since the last sum is equal to ( ^ 1 , which proves that our expression satisfies 

relation Since it also gives the expected values for r = or s = 0, this achieves 
the proof. 

Note however that those two formulas do not give any relevant information on 
the entropy since r and s the bigger, the most the tiling looks like a square lattice, 
with two defect lines (the two worms of the two other families) that only have a linear 
contribution to the entropy. Therefore the entropy per tile vanishes when r and s tend 
to infinity. 

4-2. Numerical results 

In table 0, we list some coefficients aj obtained via the method exposed above: we 
construct the oriented graph of the configuration space of the corresponding 3 — > 1 
problem and apply the generalized descent theorem, as described in the codimension- 
one case in reference Its label is attached to each (oriented) edge of the graph 
of this configuration space, as prescribed in the generalized descent theorem scheme, 
and the coefficients aj are computed recursively: if aj{v) is the number of walks in 
the configuration space from the extremal vertex O to the vertex v, with j descents, 
then the sequence (oj(f ))j=o,...,j„ax only depends on the vertices under v in the graph. 
Then the coefficients we are interested in are equal to aj{So), where Sq still denotes 
the extremal vertex associated with O. 

One checks that the corresponding sums ^ aj are in agreement with those 
computed in section 3^ (equation ^). Note that in table |l|, the symmetry aj = UM-j 



observed in codimension one is lost {M is the maximum number of descents). Indeed, 
this symmetry remains valid for each partition on tiling problem, but the maximum 
number of descents depends on the base tiling. 

The derived counting polynomials provide enumerations of tilings, as well as 
entropies per tile|f| of finite-size systems. Some examples are listed in table ||. The 
interest of the method is that it gives access to enumeration of tilings for arbitrarily 
large fc4, if fci, k2 and fca are fixed. Moreover it is technically much easier to implement 
than a brute force enumeration method, and very much faster as well, in terms of 
computational time. 

Even though it is not possible to make any reliable fit with few finite-size values, 
it is rather clear from the available data that the diagonal^ entropy converges rapidly 
to its limiting value. Note that in the 2 — > 1 case as well as in the 3^2 one, where 
exact enumeration formulas are known (see reference 1^ for a review, for instance), 
the asymptotic behavior of the finite-size corrections to the entropy can be derived : 
they decrease like log(fc)/fc. Fitting such a behavior with the numerical values, we get 

t We recall that the entropy per tile is 5 = In W/N where W is the number of configurations (tilings) 
and N the number of tiles. 

t The entropy is said to be diagonal when all the boundary side lengths fc; are equal. 
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aQ ; ai ; . . . 


2,2,2 


20 ; 220 ; 703 ; 943 ; 566 ; 166 ; 21 ; 1 


3,2,2 


50 ; 1281 ; 9775 ; 32304 ; 53175 ; 46343 ; 22095 ; 5755 ; 774 ; 47 ; 1 


2,3,2 


50 ; 1240 ; 10472 ; 40378 ; 77328 ; 75652 ; 36506 ; 7958 ; 648 ; 18 


2,3,3 


175 ; 9792 ; 183223 ; 1611390 ; 7581596 ; 20313994 ; 31942744 ; 
29678550 ; 16076840 ; 4906164 ; 794328 ; 62142 ; 2088 ; 24 


3,2,3 


175 ; 10372 ; 184113 ; 1445070 ; 5924665 ; 13826440 ; 19251677 ; 
16431348 ; 8710059 ; 2861124 ; 569191 ; 65214 ; 3943 ; 108 ; 1 


3,3,3 


980 ; 119284 ; 4736040 ; 88959048 ; 922861456 ; 5735679224 ; 22400451966 ; 
56586512056 ; 93968296600 ; 103217016568 ; 74801020694 ; 35369632364 ; 
10693166706 ; 2003702920 ; 222619576 ; 13801976 ; 439638 ; 6272 ; 32 


3,3,4 


4116 ; 990574 ; 75291817 ; 2672974232 ; 52557540678 ; 628628119744 ; 
4845859698991 ; 25007135636872 ; 88641414434386 ; 219565301033744 ; 
o841o84ooi48y98 ; 4/ /33iioo/0/ZoO ; 4214^2964232612 ; 26343i6o490o3o4 ; 
115559997005453 ; 35098071282418 ; 7238626577471 ; 987285691504 ; 
85977846450 ; 4564265102 ; 138792310 ; 2208928 ; 15936 ; 40 


3,4,4 


24696 ; 11185183 ; 1658701257 ; 117639867825 ; 4696728888239 ; 115554431503049 ; 
1855639954964533 ; 20237165017326054 ; 154261056214441072 ; . . . 


4,4,4 


232848 ; 211868010 ; 59911555328 ; 7889440518518 ; 578616346951691 ; 
26140019431942187 ; 775751817756005455 ; 15811577667366075305 ; . . . 


4,5,5 


16818516 ; 52683466776 ; 49453853710872 ; 21112489152560570 ; 
4940628646460445115 ; 704860523557345706986 ; 65676322673579106872954 ; . . . 


5,5,5 


267227532 ; 1658888888852 ; 2898208633474138 ; 2212967878070760376 ; 
903353585201401013350 ; 221402610595368245987868 ; . . . 



Table 1. Coefficients aj associated with some octagonal tiling problems. 





Number of tilings 


# tiles 


Entropy 


1,1,1,1 


8 


6 


0.34657 


2,2,2,2 


5383 


24 


0.35796 


3,3,3,3 


273976272 


54 


0.35979 


4,4,4,4 


1043065776718923 


96 


0.36022 


5,5,5,5 


296755610108278480324496 


150 


0.36031 



Table 2. Some tiling enumerations computed with the previous coefficients and 
the corresponding entropies per tile. We have only listed diagonal entropies. The 
number of rhombi is given in column 3. 

a limiting diagonal entropy close to 5 = 0.36(1). The precision of this entropy cannot 
be refined beyond the second digit with the small amount of values we have got. 

Previous entropy calculations via transfer-matrix methods where derived 
concerning octagonal tilings, but in the case of periodic or free boundary conditions 
leading to a limiting value S = 0.434. The difference between both results is due 
to the strong macroscopic effects of boundary conditions in those random tiling 
systems [|9[ p2| . 

Conclusion 

Among the large class of random tilings, this paper is devoted to fixed boundary 
codimension-two tilings of rhombi ("octagonal tilings"). We have established 
combinatorial properties of the configurational spaces of such tilings, extending results 
previously derived in the more restricted case of codimension-one tilings. Octagonal 
tilings are more closely related to real quasicrystals than are codimension-one tilings. 
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Moreover, many of the results presented here can (at least partially) be extended to 
two-dimensional tilings of rhombi of any codimension, and even to any-dimensional 
tilings. 

The present analysis provides additive formulas which simplify significantly the 
enumeration of finite-size tilings. In a geometrical viewpoint, these formulas come 
from a decomposition of the configuration space into elementary volumes, called 
normal simplices. The number of configurations in each of these simplices is known. 
But it is necessary to take into account interfaces between those volumes to avoid 
multiple counting, which is achieved by the generalized descent theorem. The number 
of simplices in this decomposition is also derived in the general two-dimensional 
case: these simplices are put in one-to-one correspondence with a class of paths in 
a configuration space of tilings of same codimension, but smaller dimension, which 
can be counted. 

The new insight on the sets of octagonal tilings provided by this analysis will be 
useful to study topics such as diffusion in these configuration spaces, which is directly 
related to the rate of convergence of fiip dynamics towards the equilibrium distribution. 
This problem has already been treated in the case of hexagonal tilings ^ |3|, Q , 
but is still an open question in higher codimension plane tilings. Significant progress 
will be published separately. It would also be of high interest to understand how the 
introduction of energetic interactions between tiles is translated in the configuration 
space and how it modifies the dynamics. Indeed, a realistic model of quasicrystals 
requires one takes into account energy, which can be in first approximation modeled by 
tile interactions; glass-like slow dynamics are likely to appear in this case |3|, even 
though no glassy behavior has been explicitely exhibited in rhombus tilings yet ||37|, ^ . 
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Appendix A. Coefficients aj and maximal vi^alks in configuration spaces 

We consider D ^ 2 fixed boundary tilings, described as generalized partitions of 
height ko on D — I 2 membranes or tilings. To each such membrane is attached 
a partition problem, and therefore a set of coefficients aj . The sum of the coefficients 
ttj of the counting polynomial of all these D ^ 2 tilings is equal to the sum, running 
over all the relevant D — 1 2 membranes, of the sums of the coefficients aj on each 
such membrane. 

Consider first a 13 — 1 —> 2 membrane, denoted by M, and a partition problem on 
Ai. According to results of section 3.1, the sum of the coefficients aj of this partition 



problem is equal to the number of walks, in the configuration space of the partition 
problem of height 1, between two extremal D ~ 1 ^ d configurations, O and Sq. All 
the parts are of O and 5*0 are respectively equal to and 1. Such a walk is denoted 
by a sequence O = Po, -Pi, ■ • ■ , Pk-i,Pk = So- 

We use again the grid representation of tilings and the formalism introduced in 
section the partition problem on A4 is seen as a partition problem on the vertices 
of the corresponding grid (which has D — 1 families), Q, called the subgrid of the 
problem. The vertices of G are ordered on each de Bruijn family. A partition Pi of 
height p = 1 consists of marking each vertex of ^ by a or a 1. The vertices marked 
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with a and those marked with a 1 are separated by the only de Bruijn line of the 
D-th family of the original grid. This latter line is denoted by So- Our goal is now 
to encode Sd by a D — 1 — > 1 grid, or in other words by a Z? — 1 1 tiling. 

Now the section of the subgrid Q by Sd is precisely of this type. Indeed, if we 
identify the intersection of a de Bruijn line of the fc^-th family and Sd with a i-tile, 
this latter section is a sequence of tiles, ki of each family i, that is to say a 13 — 1 ^ 1 
tiling. 

Therefore the sequence (Pi) is coded by a sequence {d) of such tilings, as 



illustrated in figure |A1 
g of this D 



Since the same argument can be applied to any subgrid 
d problem, all the walks counted by the sum of the coefficients Uj of 
the counting polynomial can be seen as such sequences (Ci). 





Figure Al. Left: A _D — 1 — > 2 grid and the successive sections d by a line 
of the D-th family. Each section can be seen as a D — 1 — > 1 tiling where the 
different species of tiles are represented by dots of different colors. Two successive 
sections differ by a single flip, that is the exchange of two tiles. In this example, 
D = 4, fci = 3, fc2 = 1 and = 2. Right: The corresponding labeled dual tiling. 



It is now rather clear that all the configurations Cq — O and Ck — Sq are the 
same for the partition problems on all subgrids Q: they are the two configurations 
where all the tiles of a same family are adjacent and where the families are ordered 
according to the grid configuration at infinity, which does not depend on Q, but only 
on the vectors it^. These two configurations will be denoted by Cmhi f^nd Cmaxi as in 
figure Al. Likewise, two successive tilings of a sequence, Ci and Ci+i, only differ by 
a tile flip, that is the exchange of two adjacent tiles, since the two successive sections 
"surround" a single vertex of Q. And there is a natural order between two successive 
such tilings: considering how tiles are ordered on Cmin and Cmax, it is clear which 
tiling should come first in the sequence. 

Conversely, let us establish that such a sequence (Ci) of 13 — 1 ^ 1 tilings, going 
from Cmin to Cmax, and where two successive tilings only differ by a single flip and 
respect the above order, contributes towards the sum of the coefficients aj of a unique 
partition problem on a 13 — 1 ^ 2 membrane - and therefore towards the sum of the 
coefficients aj of the global 13 ^ 2 problem. The proof is rather straightforward and 
is also illustrated in figure Al: considering two successive tilings, Ci and C^+i, the 



two dots that represent the two flipping tiles of different families are joined by two 
crossing segments; the so-obtained vertex is labelled by i -I- 1; all the other tiles are 
joined by horizontal segments. Then one reconstructs a complete 13 — 1 ^ 2 grid, the 
vertices of which are labeled by numbers increasing on each de Bruijn line. This is 
precisely the kind of object counted by coefficients aj. 

Note that all the walks counted by ^ aj have the same length, since all the 
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subgrids Q have the same number of intersections. 

Appendix B. Number of mziximal walks in a n ^ 1 configuration space 

In this section, we derive the number yl„(fci, . . . , A:„) of walks {Pi)i=Q^...^K in the n — s- 1 
configuration space. Note that for sake of simphcity, we note n instead of 13 — 1. We 
shall prove that: 

l<i<j<n 



n 



where Kij = 2(fci+i + . . . + kj-i) when i < j. 

We giv e two proofs: a purely al geb raic one in the general case, using results by 
Stanley , and a combinatorial one |2^ . Note that the case where all the parameters 
ki are equal to 1 was already treated by Stanle y |39|] , and that Edelman and Greene 
derived a nice combinatorial proof in this case |t40| . We suppose that the notions of 
Young tableaux and standard Young tableaux are known^ 

But before all, we need to introduce the relation between the tilings considered 
in this paper and a class of computer science objects, the so-called sorting algorithms. 
This analogy will help the presentation of the algebraic proof and will be useful in the 
combinatorial proof. 

Appendix B.l. Tilings and primitive sorting algorithms 

In the sorting language, a comparator [i; j] acts on a list [xi, X2, ■ ■ ■ ^ Xn) of numbers as 
follows: Xi and Xj are respectively replaced by min(a;i, Xj) and max(a;i, Xj). Following 
Knuth jH], we call a complete sorting algorithm a sequence of such comparators 
which sorts in the increasing order any list of real numbers (xi, X2, ■ ■ ■ , Xn)- This 
sorting algorithm will be called primitive if each comparator can be written [i, i + I]. 
We also suppose that this algorithm is not redundant, that is to say it does not 
contain any comparator [i,j] that could be suppressed because previous comparators 
already insure that Xi < Xj . Knuth shows that a sequence of comparators is a sorting 
algorithm if it correctly sorts the completely reversed list (n, n — 1, . . . , 1). This means 
that a complete primitive sorting algorithm is a sequence of comparators [i, i + 1] that 
transforms the list (n, n — 1, . . . , 1) into the list (1,2,..., n). 

Such an algorithm can have a diagrammatical representation as follows: the 
n variables Xi are represented by n horizontal lines. Each comparator + 1] is 

represented by a crossing ^XZ between lines i and i + Figure Bl illustrates 



this construction. A continuous line follows a number during the sorting process. For 
example, the greatest number is on the top at the beginning and in the bottom at the 
end. Since every number must be compared to every other one, and since there is no 
redundancy, there are (2) crossings. 

f A Young tableau is a stacking of square boxes on a line, as displayed in figure |B4| (left), with the 
only constraint that the number of boxes in columns decreases from left to right. The shape of the 
tableau is the decreasing sequence of the column heights. A standard Young tableau of a given shape 
is simply a numbering of the cells of the tableau by integral numbers, running from 1 to the number 
of cells and increasing in rows and columns. Figure |B4| provides examples. For a presentation in 
relation with the symmetric groups Sn, the reader can refer to Jill. 
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Figure Bl. A diagram associated with a sorting algorithm acting on five element 
lists. A line follows a number during the sorting process. Each pair of lines cross, 
only once. The two circled comparators can be exchanged without changing the 
corresponding tiling. 

We are now able to establish the link between those algorithms and n ^ 2 tilings, 
more precisely with their de Bruijn representation. Indeed, the analogy between the 



diagram of figure Bl and a de Bruijn grid with one line per family is straightforward: 
each continuous line of the diagram represents a de Bruijn line and crosses exactly 
once every other line. 

However, there is a fundamental difference between both systems: different 
sorting algorithms can represent the same de Bruijn grid since only the crossing 
topology is meaningful. For example, in figure ^l|, the fourth and the fifth comparator 
(i.e. [4,5] and [2,3]) are applied in this order (these comparator are circled in the 
figure). If they were applied in the reverse ordre, the algorithm would be different 
whereas the de Bruijn grid would be the same. 

Therefore we are led to define equivalence classes of sorting algorithms [^, |2^ . We 
say that two comparators [i, z + 1] and [j, j + 1] commute if |i — j| > 1. Two algorithms 
are equivalent if they differ by a finite number a comparator commutations. These 
equivalence classes of n-element sorting algorithms are in one-to-one correspondence 
with n-family grids with one line per family, and therefore with tilings inscribed in 
polygons of side 1. The number An of equivalent classes has been computed by 



Stanley [Q (see also Edelman and Greene jS^]) and is given by equation B.l in the 
case where ki = 1 for all i: 

-^"(^'^'■■■'^^= l»-i3»-ffU-3)^ - ^^-'^ 
We need to generalize this point of view to systems with more than one line per 
de Bruijn family, which leads to the definition of partial sorting algorithms. These 
algorithms are related to pre-sorted lists of numbers. Indeed, let us suppose that 
we have n families of ki numbers each (i = l,...,n), and that in each family, the 
numbers are already pre-sorted in the increasing order. Then we are interested in 
the algorithms which order the whole set of these numbers in the increasing order. 
We call them partial sorting algorithms. The ideas are essentially the same as in the 
previous case, except that, since the numbers of a given family are already ordered, the 
corresponding lines do not need to cross. The corresponding diagram is similar to a de 
Bruijn grid with n families of lines, hi lines in each family. The tilings are equivalence 
classes of such algorithms. They are inscribed in polygons of sides /ci, . . . , A:„. 

In the partial sorting case, the reference list to be reversed is not (n, n — 1, . . . , 1) 
any longer but a list wq were some elements are already sorted. If = ki + k2 + ■ ■ ■ + ki, 
then 

Wo = (Hn-l + ^, ■ ■ ■ , Kn-l + kn, (B.3) 
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,1,2,, 



,fci). 



There are n pre-sorted blocks of ki elements each. An example is provided in figure B2 




Figure B2. A partial sorting algorithm in the case of n = 3 families, containing 
respectively fci = 3, fc2 = 1 and = 2 pre-sorted elements. The sequence uiq 
appears vertically on the left of the figure. 



Note also that a sorting algorithm can be seen as an ordering (or a labeling) of 
the comparators of the corresponding equivalence class, since the comparators come in 
a natural order: the labels arc increasing on each de Bruijn line, from left to right. As 
a consequence, sorting algo rith m s are in one- to-one correspondence with labeling of 
tilings, as d efin ed i n section 3^ or Appendix A . This point is striking when comparing 
figures Al and B2. Therefore the sums of the coefficients of a D — > 2 problem is 
equal to the number of partial sorting algorithms of the suitable pre-sorted D — 1 
families of variables. This point will be helpful in the following section. 



Appendix B.2. Walks and symmetric groups Sn (algebraic proof) 

The main result used in this section is a theorem by Stanley (theorem 4.1 of 
reference and its corollary 4.2). We first need to expose these results and to 
translate them in terms of our notations and definitions. It is the object of the first 
paragraph. Then, using the above equivalence between coefficients Oj and partial sorts 



we shall apply Stanley's theorem to the case of interest here, and derive relation B.l 



Appendix B.2.1. Stanley's theorem: The symmetric group Sn of permutations on n 
elements is generated by the transpositions cji = {i^i -\- \): any permutation w can be 
decomposed in products of transpositions. A decomposition is said to be minimal if, 
using the relations between the generators ai: 

crf^l and a.,ai+ia.i = cr.,+iCJiai+i, (B.4) 

it cannot be simplified into a shortest decomposition; then all the reduced 
decomposition of w have the same length||], denoted by l{w). w = 

i ^ ^ y we define 

\ ai 02 ■ ■ ■ On J 

ri{w) — Card{j : j < i and oj > flj} (B.5) 
Si{w) — Card{j : j > i and oj < ai}. (B.6) 

f The length of a decomposition is the number of generators necessary to write it. 
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Then X{w) is the sequence obtained by arranging the numbers ri{w) in descending 
order (and ignoring any O's); fJ.{w) is the conjugate^ to the sequence obtained by 
arranging the numbers Si(w) in descending order. 

Stanley ]39| states that: if X{w) = fJ'i'w), then the number of reduced 
decompositions of w is equal to the number f^^^^ of standard Young tableaux of shape 
\{w). 

Appendix B.2.2. Application to walks inn 1 configuration spaces: We still denote 
by Wq the permutation as defined in the previous section. The idea is to see each 
comparator [z, z + 1] as a generator Ui = (z, z + 1) of a reduced decomposition and to 
identify a sorting algorithm with a reduced decomposition of wo- More precisely, the 
set of comparators can be seen as the gen erat ors of a group, and they obey the same 
relations as the Ci: they satisfy relations B^, and [z,z + 1] and [j, j + 1] commute if 
|z — j| > 1. Finally, the non-redundant character of sorts is equivalent to the reduced 
character of decompositions. As a consequence, the number A„(fci, . . . , of partial 
sorts is equal to the number of reduced decompositions of wq- 

Now, in order to apply Stanley's statement, we must compute the quantities 
\[wq) and /i(wo), as defined in reference |^^: in the present case, one gets 

n(wo) ~ kn + ■ ■ - ^-kj if kn + ■ ■ ■ + kj + l>i>kn + . . . + fcj-i(B.7) 
St{wo) — if kn + ■ ■ ■ + kj + I > i > kn + ■ ■ ■ + kj-i. (B.8) 

X{wo) = ^J.iwo) = {kn + . . . -I- /C2, fcn + ■ ■ • + fcs, ■ ■ ■ 



Thus 



ki times 



times 



k^ ), (B.9) 

/Cti-i times 



and Stanley's theorem applies. The Young tableau of shape A(wo) is represented in 
figure B3. Such a tableau will be called a block tableau of size (fci, k2, ■ ■ ■ ,kn) in the 
following. The number v4„(/ci, k2, ■ ■ ■ , kn) of reduced decompositions of wq is equal 
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Figure B3. A Young tableau of shape A(uio) = /i(«)o)- The numbers ki denote 
the number of rows and columns in each rectangular block. Such a tableau is 
called a block tableau of size (fci, k2, ■ . . , k„). 



f If a decreasing sequence is drawn like a Young tableau, its conjugate is the decreasing sequence of 
row lengths. For example, the conjugate to (3,2,2,1) is (4,3,1). 
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to the number f^(''"o) gf standard Young tableaux of shape X{wq). This number can 



be derived from Young's hook-length formula |43, given a shape A, the hook 
associated with a given cell c of the tableau is the set of cells above and at the right 
of c, including c itself (figure B4). It is denoted by He- The hook length he is the 
number of cells in He- Then the number of standard Young tableaux of shape A 
is: 



where N is the total number of cells and the product runs over all cells. 
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Figure B4. Left: a four-column Young tableau of shape A = (5, 3, 1, 1). Middle: 
an associated standard tableau. The hook corresponding to the circled entry is 
represented. Its length is he = 5. The hook rank of this cell is Tc = 1 and its hook 
height is 6c = 4: this tableau is not balanced. Right: a balanced stair tableau of 
order 5. 

This hook-length formula can know be applied to the above tableau: the 
hook- length of the cell in the upper right corner of block ki x kj is Kij + 1, 
from which one deduces all the hook-lengths on the block and their product: 
(K- +k- + k- ~ IV^'^UK-- - 

{K',+k.-r)miK,,+k,-l)m - ^ - T..<,hk,, from which equa- 



tion B.l follows 



Appendix B. 3. Combinatorial proof 

In this section, we provide a combinatorial bijective proof of the previous result pl| 



This proof follows the same scheme as the proof by Edelman and Greene [40 in the 
case where all the parameters ki are equal to 1. We need first to introduce the notion 
of balanced tableaux 

We consider a Young tableau, together with a labeling of its cells, running from 
1 to the number of cells, K. Note that now this tableau is not necessarily standard, 
that is to say the labels are not necessarily ordered in each row and in each column. 
Given a cell c and its hook He, we define the hook rank rc of c as the number of cells 
of He whose labels are smaller or equal to the label of c. We also define the hook 
height 9i as the number of cells above c (including c itself) . The tableau is said to be 



balanced if for all cells c, rc = 0c (see figure B4). 

At last, a tableau of shape (n — 1, n — 2, . . . , 1) (figure B4, right) will be called 
a stair tableau (of order n). Note that a stair tableau is a particular case of block 
tableau, as defined in the previous section. 
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In the case fcj = 1, the situation is as follows: Edelman and Greene build a 
bijection between complete sorting algorithms on n elements and balanced block 
tableaux of order n and then a bijection between those balanced tableaux and standard 
tableaux of order n, which are then enumerated via the hook-length formula. In the 
following, we use and generalize these results to the case of partial sorts. In fact, 
in reference the authors also establish the bijection between balanced tableaux 
of any shape and standard tableaux of the same shape. As a consequence, we only 
need to generalize the first bijection between partial sorting algorithms on n families 
containing ki elements each and balanced block tableaux of size (fci, . . . , /c„), as in 
figure Once this correspondence is established, the rest of the proof is based upon 
the hook-length formula, as in the previous algebraic proof. The rest of the section is 
devoted to this correspondence. 

To begin with, let us consider a complete sorting algorithm on n lines: it is a 
sequence of crossings between de Bruijn lines. These crossings are labeled by integers 
running from 1 to K = (2)1 from left to right. Following Edelman and Greene, the 
intersection label of two lines indexed by a and b is denoted by tab- In the stair 
tableau, this number is written in the cell situated on the a-th column and on the 
{n — b + l)-th line (starting from the bottom): in figure B4 (right), the so-obtained 



tableau corresponding to the complete sorting algorithm of figure Bl is represented. 
Then it can be proven that this tableau is balanced and more precisely that this 
construction establishes a bijection between both classes of objects. 

Let us now focus on partial sorting algorithms. Since there are couples of de Bruijn 
lines which do not intersect, it is rather natural to consider block tableaux, which are 
stair tableaux where some cells are missing. More precisely, we will consider partial 
sorting algorithms and block tableaux as amputated complete sorting algorithms and 
amputated stair tableaux, respectively. The idea is to define canonical amputations 
in order to preserve Edelman and Greene's bijection between amputated objects, as 
discussed below. 

Figure p35| illustrates the tableau amputation process: given n integers /ci, . . . , fc„ 
and a stair tableau of order N — ^ki, n small stair tableaux of order ki are removed 
from the large one in order to get a block tableau of size (A^i, ■ ■ ■ 5 ^n)- 



Figure B5. Amputation of a stair tableau. 
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As far as the amputation of sorting algorithms is concerned, we first need to 
define canonical complete sorts: the simplest way to characterize them is by their 
corresponding stair tableaux, the cells of which are increasing from left to right and 
from bottom to top, as illustrated in figure p36| . They are usually referred as "bubble- 
sorts" in the literature. 
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Figure B6. 

tableau. 



A canonical complete sort and the corresponding canonical stair 



If we consider now a partial sorting algorithm with ki elements in each family, 
it can be canonically transformed into a complete sort: we simply add n canonical 



complete sorts on ki elements at its end, as in figure B7. These sorts appear in the order 
of their indices i. The so-obtained sort is very particular since it ends with n canonical 
sorts. Therefore it will be called a sequential complete sort of order (fci, . . . , A:„). By 
construction, there is a one-to-one correspondence between sequential complete sorts 
of order (fci, . . . , and partial sorts of the same size. 
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Figure B7. A complete sequential sorting algorithm: it is a partial sorting 
algorithm followed by n canonical complete sorts on ki elements. 

Let us now characterize the structure of the Young tableau t associated with 
such a sequential sort. The (^^2 ^ + ('2^) ■ ■ • ^ ('2") intersections in the sort 
are those appearing in the n canonical complete sorts. It is easily checked that the 
corresponding labels in the large stair tableau appear in the n small canonical stair 
tableaux involved in the amputation process. As a consequence, the labels of the 
amputated tableau run from 1 to K, where K is its number of cells. These labels code 
the K intersections of the partial sort remaining of the original sequential complete 
sorting algorithm. 



To sum up, as it is illustrated in figure |B8| , starting from partial sorts, we 
biunivocally construct complete sequential sorts, then balanced tableaux, the K last 
labels of which are situated in the n stair sub-tableaux. When these sub-tableaux are 
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removed from the large one, we obtain a class of block tableaux of size (fci, . . . , A:„), 
which will be called pre-balanced tableaux. Remember now that are goal is to establish 
a bijection between partial sorts and balanced block tableaux. Thus we need to 
construct a bijection (denoted by R in figure BS) between those pre-balanced tableaux 
of size (fci, . . . ,kn) and balanced block tableaux of the same size. 



Balanced block tableau 



Complete sequential sort 



Balancetl stair tableau 



Figure B8. Sequence of bijections establishing tlie one-to-one correspondence 
between partial sorting algorithms on n families with ki numbers each and 
balanced block tableaux of size (fci , . . . , In the lower right stair tableau, 

the grayed sub-tableau is a pre-balanced block one. The bijection R puts such 
tableaux in one-to-one correspondence with balanced block tableaux of the same 
shape. 



Bijection R: actually, R is an involution: in each group of ki columns, it inverses 
the order of columns; for example, the first column becomes the fci-th one, the second 
one becomes the (fci — l)-th one and so on. Likewise, the (fci -I- l)-th column becomes 
the (^2 — l)-th one. There is a rather deep reason for such a column permutation: 
the n canonical complete sorts added at the end of a partial sort in order to make 
it complete also reverse the order of lines in each group of ki lines. The role of R is 
to keep track of this fact. It is now a rather technical task to prove that R provides 
balanced tableaux and that it is a bijective map. 

We shall temporarily admit the following results which will be proven at the end of this 
section: in a pre-balanced tableau, in each block ki x kj, the labels are decreasing in each 
line and column (from bottom to top); in a balanced block tableau, they are increasing in 
lines and decreasing in columns in such a block. 

Consider a pre-balanced block tableau and in this tableau a hook associated with 
the cell c situated in the block ki x kj, and in this block in the line u (from bottom to top) 
and column v {1 < u < ki and 1 < v < kj). This hook comes from a larger one, H'^, in the 
stair tableau which has been amputated. In the process the hook has lost C cells, Ci on its 



right and C2 on its top (see figure B9). Thus the height of He is equal to 



9c = O'e — C2- 
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On the other hand, all the lost cells, coming from the n removed stair tableaux, had labels 
larger than the label of the cell c. As a consequence, the hook rank remains unchanged: 

Tc = r'c. 




Figure B9. A hook He of a prc-balanccd tableau. As compared to the original 
hook H'^ (before the amputation), it has lost cells on its top and on its right 
(dashed) . By the involution R, c receives the label of the cell c'^ . 



Now, R permutes the labels of the block tableau but does not change its shape. In the 
process, the cell c receives the label of the cell c° still situated in the line u, but in the column 

kj — V + 1. 

For a given quantity A assigned to each cell, we denote by A" the value of this quantity 
for the cell c" in the pre-balanced tableau, and its value in its image by R. In particular, 
= 9", and wc have just proven that r" - 9" = . Moreover, rf = r" - (kj - v°). Indeed, 
by R, the hook He receives the labels of the hook Hcq outside the block ki x kj; and in 
the block, the labels, which were decreasing in the line u of the pre-balanced tableau, are 
increasing in its image by R. As a consequence, 

R qR ^0 /, On 

rc. - = ^2 - [kj -V )■ 

Now C2 ~ {kj —v"^) by definition. Therefore — 9^ and the tableau is balanced. Conversely, 
given a balanced tableau, since the labels are increasing in each line of each block, one proves 
that its image by the involution R is pre-balanced. We have established the bijection. 

We need to prove the two above assertions about the order of labels in blocks of balanced 
and pre-balanced tableaux. We only give sketches of the proofs. As far as pre-balanced 
tableaux are concerned, the proof is rather straightforward: a block ki x kj contains labels 
associated with all the intersections of two families of lines. If those lines are isolated from the 
rest of the tiling, it becomes clear that the order in which intersections occur is constrained. 
For balanced tableau, the proof is more complex. The basic idea is to construct a proof by 
"planar induction" : we prove that if a suitable V property is true for cells above and at the 
right of a given cell c, then it is also true for c. Then if V is true for cells on the upper right 
corners of the tableau, it will be true for every cell. In the present case, if tc still denotes the 
label of the cell c, then the property reads: 

^(c); in the block ki x kj to which c belongs, the cells above c have labels smaller than tc, 
while the cells at its right have labels greater than tc ■ 
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A proper use of the balanced character of a block tableau proves that it satisfies the above 
planar induction principle. 

Appendix C. Proofs of section ^ 

In this appendix, we only give sketches of proofs for the results used in section ^. The 
complete proofs can be found in the relevant references. 



Appendix C.l. Elnitsky's formulas 

In order to prove relations ^ and 0, we need, in the octagonal slightly different 

representation of de Bruijn grids than those presented in the introductory section: as 
displayed in figure the lines of the two first families form a square lattice (of sides 
r and s in this case), on which the lines of the third and fourth de Bruijn families 
run: they are represented by directed walks on the lattice, going from one corner to 
the diagonally opposite one. According to whether the two sides of length one are 
adjacent or not, the two paths have the same starting and ending points (right) or not 
(left). In this representation, the de Bruijn line intersections must be distinguished 
to avoid possible ambiguities due to path tangency (it will be called the distinguished 
vertex in the following). For e xam ple, the octagonal tilings of figure 11 are represented 
by the pair of paths of figure jCl]. 
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Figure CI. Two lattice-path diagrams, associated with the two tilings of 
figure hll. The de Bruijn fines of the two single-line families are represented by 
a pair of paths running on a square lattice. These paths go from one corner to 
the diagonally opposite one. In case of ambiguity, the vertex where the de Bruijn 
lines intersect is distinguished by a circle. There are two cases: either the two 
sides of length one of the polygonal boundary are adjacent (right) or not (left). 

The proof of relation |^ is now straightforward: if (a, b) are the coordinates of 
the distinguished vertex in the grid, then the number of configurations is the product 
of the 4 number of choices for the four pieces of paths going from (a, b) to the four 
corners, that is the product of four binomial coefficients. Now summing over all (a, b) 
configurations one gets relation ^ 

The proof of relation ^ is a little more complex and involves some modifications 
of the pair of paths, as usual in this kind of calculation. The idea is to exchange the 
two paths after the distinguished vertex, in order to get non-crossing (but possibly 
touching) paths, and then to shift some parts of those paths in order to get a non- 
touching pair. The latter pairs can be counted with help of the determinental Gessel- 
Viennot method The interested reader will refer to reference for more details. 
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Appendix C.2. Brock's recursion relation 

Let us now focus on the recursion relation We use again the lattice-path 

representation, as defined in the previous section. Given a pair p of crossing paths, 
with no a priori distinguished vertex, let us denote by Ip the length of their intersection 



(see figure C2). Note that this intersection can be a point [Ip = 0), a horizontal or a 



vertical segment (of length Ip). If Vr.s is the set of pairs, then 



r 

since there are Zp + 1 possible choices for the distinguished vertex and there are 
/ r + s V 

( ^ 1 such pairs of paths. 
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Figure C2. When the two path intersection is a segment (vertical or horizontal, 
of length Ip > 0; here Ip = 3), the same pair can represent several tilings, 
depending on the position of the distinguished vertex. If this point is not the 
last possible one (i.e. the point with the highest coordinates), then a strip of the 
lattice can be removed (dashed in this figure, at the right of the black distinguished 
vertex). The so-obtained diagram is associated with a smaller tiling (r — 1 X s in 
the present case; rxs — lifthe segment were vertical) . 



In the case where Ip > 0, if the distinguished vertex is not the point with the 



highest coordinates {Ip possibilities), as in figure C2, then a vertical or horizontal strip 
of the square lattice can be removed without changing the nature of the diagram: it is 
the strip of width 1 , at the right of (respectively above) the distinguished vertex if the 
intersection is horizontal (respectively vertical). As a result, one still have a square 
lattice (but one of its side lengths is lowered by 1) with a pair of paths (but the length 
of their intersection has been lowered by one in the process) . That is why 

E + E (^p + i)- (c-2) 



In conclusion, 
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r + s 
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- E ('p + i)+ E (^p + i) 
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r + s 
r 



which achieves the proof. 

Appendix D. What about the general D ^ d case? 

In this last appendix, we discuss what we know and what we do not know about the 



general D ^ d case which was more widely studied in reference 21 . The different 
points and results tackled in the present paper are discussed. It will become clear 
that the possible existence of cycles in the partition-on-tiling problems is the major 
obstacle to simple generalizations. 

Appendix D.l. Partition- on-tiling point of view and configuration space 

As far as partitions on tilings are concerned, all that as been said in the octagonal case 
can be transposed to the general case: a fixed boundary D ^ d tiling can be 

coded in a single way as a partition on a D — 1 — > d tiling. The de Bruijn lines are still 
defined as lines joining together the middles of opposite faces of rhombic tiles, and the 
parts of those partitions are still increasing along such lines. The reader can refer to 
Bailey |^ for a more formal treatment of this question. Note that in this case, there 
also exist de Bruijn families of hyper-surfaces, associated with an edge orientation. 

However, as it was suggested in reference [||, the geometry of configuration spaces 
might be more complex beyond the octagonal case. Indeed, among the order relations 
Xi > Xj between the parts of a partition-on-tiling problem, nothing forbids a priori 
the existence of cycles of inequalities, such as > Xi^ > . • . > Xi^ > Xi-^ , which 
enforces all these variables to be equal. 

At least two examples of 6 — > 3 tilings (in relation with 7^3 tiling problems) are 
known which display such cycles. The first one can be found in reference [ p6| (example 
10.4.1) and the second one in (example 3.5)||. In these examples, the tilings are 
defined by their dual de Bruijn grids - which are families of 2-dimensional de Bruijn 
surfaces in a 3-dimensional space - together with an orientation of de Bruijn lines. 

When such a cycle exists, all the parts of the cycle have a "collective behavior" , 
which is not compatible with the previous description: they behave like a single 
effective part. In particular, the number of effective parts in this partition problem, 
denoted as K' , is strictly smaller than K. Thus the counting polynomial of this 
partition problem becomes: 

j=0 ^ ^ 

The counting polynomial of the whole tiling problem is a sum of such polynomials, 
with possibly many different K' . 

Moreover, the existence of cycles invalidates the proof of the connectivity of the 



configuration space (section |2^) . As far as we know, this point is an open question 
in the general D d case. Note however that whenever one can prove that order 
relations on fibers contain no cycles, then the configuration space is connected. 

f Note that in these references (proposition 10.5.7 of [ p^ and corollary 4.5 of ji?)), it is also stated 
that no such cycles exist in two-dimensional tilings. 
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Appendix D.2. Decomposition in simplices - descent theorem 

We prove that the existence of cycles does not alter the previous results about the sums 
of coefficients aj and walks in configuration spaces, provided these objects are suitably 
defined: if we focus only on coefficients af associated with configuration spaces of 
effective dimension K (and not K' < K), then the sum of these coefficients is equal 
to the number of maximal walks in the corresponding Z) — 1 — > d — 1 configuration 
space. Note that this quantity still characterizes the leading coefficient of the 
counting polynomial (of degree K) as p goes to infinity. 

More precisely, if there exists a cycle in the partition problem on a membrane, we 
have just seen that the coefficients of this problem do not contribute to S^. On 
the other hand, let us consider a step in a walk in the 13 — 1 — > d— 1 configuration space: 
we have seen that the part x which differentiates the two consecutive configurations 
of this step is such that all parts lesser (resp. greater) than x in the graph are equal 
to 1 (resp. 0). This part is equal to in a partition and to 1 in the other one. As 
far as the parts of the cycle are concerned, since they are all equal, they cannot but 
jump from to 1 all together, which is not conform to our definition of walk in the 
configuration space, in terms of single elementary flips. Conversely and for similar 
reasons, a maximal walk in the configuration space cannot give a membrane with 
cycles. 

In conclusion, as well as the number of walks are not concerned by partitions 
with cycles, and the above result remains valid. 

Before going on, let us specify what the extremal tilings become in larger 
dimension: Cmin and Cmax are defined by partitions where all the parts are equal to 
and 1, respectively. Therefore the corresponding tilings present a faceted aspect, as 



on figure Dl, Note that among the different possible faceted tilings, the two extremal 



ones depend on how the D-th de Bruijn family of surfaces is chosen. 




Figure Dl. Two examples of extremal 4 — > 2 tilings displaying a macroscopic 
faceting. Such tilings are extremal configurations for walks associated with a 
5 — > 3 problem. There are four similar pairs of faceted tilings which could be 
chosen as extremal configurations, depending on how the fourth de Bruijn family 
is chosen among the four possible ones. 

At last, we recall that there cannot exist a descent theorem as simple as the 
octagonal one in the general case. Indeed, its derivation is closely related to the 
existence of a zero-descent simplex in each partition-on-tiling problem, that is to the 
if-dimensional character of the associated configuration space. But we have shown 
that this point is not granted in general, since there can exist partition-on-tiling 
problems for which the configuration space has a dimension K' smaller than K. 
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